Dissipative dynamics of a vortex state in a trapped Bose-condensed gas. 
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We discuss dissipative dynamics of a vortex state in a trapped Bose-condensed gas at finite tem- 
perature and draw a scenario of decay of this state in a static trap. The interaction of the vortex 
with the thermal cloud transfers energy from the vortex to the cloud and induces the motion of the 
vortex core to the border of the condensate. Once the vortex reaches the border, it immediately 
decays through the creation of excitations. We calculate the characteristic life-time of a vortex state 
and address the question of how the dissipative dynamics of vortices can be studied experimentally. 
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The recent successful experiments on Bose-Einstein 
condensation (BEG) in trapped clouds of alkali atoms 
|Q-|[ have stimulated a great interest in the field of ultra- 
cold gases 0] . One of the goals of ongoing studies is to 
investigate the nature of a superfluid phase transition 
in ultra-cold gases and to make a link to more compli- 
cated quantum systems, such as superfluid helium. Of 
particular interest is the relation between Bose-Einstein 
condensation and superfluidity. However, being the most 
spectacular manifestation of the phase transition in '*He, 
superfluidity has not yet been observed in trapped gases. 
A promising way of studying superfluidity in trapped 
gases is the creation of quantum vortices, as quantization 
of circulation and the related phenomenon of persistent 
currents are the most striking properties of superfluids. 

A widely discussed option of creating vortices in 
trapped gases assumes the rotation of a slight asymmetry 
of a cylindrical trap after achieving BEG, or cooling down 
the gas sample below the Bose-condensation temperature 
in an already rotating trap y,|| . Another possibility is a 
rapid quench of a gas sample near the critical tempera- 
ture, which should lead to creation of vortices even in a 
non-rotating trap H . It is worth mentioning the ideas to 
create the vortex state in a Bose-condcnscd gas by opti- 
cal means [g|,|[ , and the idea to form vortex rings in the 
regime of developed turbulence flQ]. The spatial size of 
the vortex core in the Thomas-Fermi regime is too small 
to be observed, and for visualizing the vortex state it is 
suggested to switch off the trap and let the cloud ballis- 
tically expand. Then the size of the vortex core will be 
magnified approximately by the same factor as the size 
of the expanding condensate lilt] . 

Similarly to the recently studied kink-wise condensates 
P, p^p^ , vortices are the examples of macroscopically 
excited Bose-condensed states. In a non-rotating trap the 
vortex state has a higher energy than the ground-state 
Bose condensate, i.e. the vortex is thermodynamically 
unstable ||l4|-[16|]. On the other hand, a quantum vortex 
with the lowest possible circulation (the vortex "charge" 
equal to 1), is dynamically stable (small perturbations 
do not develop exponentially with time; see |l^,fl| and 
refs. therein). Therefore, the vortex state can only decay 



in the presence of dissipative processes. 

In this Letter we discuss dissipative dynamics of a vor- 
tex state in a trapped Bose-condensed gas at finite tem- 
peratures and show how the interaction of the vortex 
with the thermal cloud leads to decay of the vortex state 
in a static trap. According to our scenario, the scatter- 
ing of thermal excitations by a vortex provides the energy 
transfer from the vortex to the thermal cloud and induces 
the motion of the vortex core to the border of the con- 
densate, where the vortex decays by creating elementary 
excitations. We calculate the characteristic life-time of 
the vortex state are discuss how the dissipative dynam- 
ics of vortices can be studied experimentally. 

We first briefly outline the main features of the vortex 
behavior in a superfluid, known from the studies of liquid 
helium. The motion of a vortex in a superfluid of density 
Ps satisfies the Magnus law (see p8[-p0{ and refs therein): 



Psi^L - V5) X K ■ 



(1) 



Here v^ is the velocity of the vortex line, and V5 the 
velocity of superfluid at the vortex line. The vector k 
is parallel to the vortex line and is equal to the circu- 
lation carried by the vortex. The force F acting on the 
vortex originates from the mutual friction between the 
normal component and the moving vortex line, and is 
usually small. Assuming the absence of friction (F = 0), 
the vortex moves together with the superfluid component 
(vl = V5). The superfluid velocity vs(r) in the presence 
of a vortex at the point Fq satisfies the equations 



rotvg = 27rK(5(r — Tq); divvs = 



(2) 



and is related to the phase of the condensate wave- 
function as vs = V0. This leads to quantization of the 
circulation: k — Zh/m [pi| , where Z is an integer (the 
charge of the vortex) and m is the mass of the condensate 
particle. Below we will consider vortex states with Z — 1, 
which are dynamically stable ( |l^,Ol and refs. therein). 
Eqs.(||) are similar to the equations of the magneto- 
static problem, with the magnetic field replaced by the 
velocity vg and the electric current replaced by k. The 



velocity field around an infinitely long straight vortex line 
is analogous to the magnetic field of a straight current: 



V5(r) = [kx r]/r 



(3) 



The vortex itself can experience small oscillations of its 
filament, characterized by the dispersion law |w(fc)| — 
Kk^ln{l/ak)/2 [^, where k is the wave vector of the 
oscillations, and a the radius of the vortex core. In a 
weakly interacting Bose-condensed gas the core radius is 
of order the healing length a = {fi /m^)^''^, where fi is 
the chemical potential. 

We will see that the dissipative dynamics of a vortex 
state is insensitive to the details of the density distri- 
bution in a gas. The spatial size of the Thomas-Fermi 
condensate trapped in a harmonic potential of frequency 
u; is R = {2fj,/mijj^y^'^. Therefore, for finding the super- 
fluid velocity vg in this case, we may consider a vortex 
in a spatially homogeneous condensate in a cylindrical 
vessel of radius R, with the vortex line parallel to the 
axis of the cylinder. For the vortex line at distance xq 
from the axis, the velocity field can be found by using 
the "reflection" method |Q. In a non-rotating trap, in 
order to compensate the normal component of the veloc- 
ity field (pi) everywhere on the surface of the cylinder, we 
introduce a fictious vortex with opposite circulation on 
the other side of the vessel wall, i.e. at distance R'^/xa 
from the cylinder axis. At the position of the vortex the 
"reflection" induces the velocity 



vs 



[l^ X Xq] 

R^-xl 



(4) 



As vs ~ vl, the vortex line will slowly drift around the 
axis of the trap. A characteristic time responsible for the 
formation of the velocity field (Q) is tr '^ R/cg, where 
Cg = \/ [ijin is the velocity of sound. Sufficiently far from 
the border of the Thomas-Fermi condensate, i.e. outside 
the spatial region where R — xq <S^ R, the drift period is 
Tdr ~ xq/vs ~ R^/k and greatly exceeds the time t/j : 

tr R a ^ 

— p « 1- 

This means that we can neglect the retardation effects 
and, in particular, the emission of phonons by the mov- 
ing vortex. In other words, the "cyclotron" radiation is 
prohibited, since the wavelength CsTdr of sound which 
would be emitted exceeds the size R of the condensate. 
According to the above mentioned magnetostatic anal- 
ogy, in a non-rotating trap the potential energy of the 
system (vortex plus its reflection) can be thought as the 
energy of two counter flowing currents. Since the currents 
attract each other, the energy is negative and decreases 
with displacing the vortex core towards the wall. In other 
words, it is energetically favorable for the vortex to move 
to the border of the vessel. Near the border the velocity 
of the vortex exceeds the Landau critical velocity, and in 
a homogeneous superfluid the vortex decays through the 



creation of phonons [g8[. In a trapped gas the conden- 
sate density strongly decreases near the border, and the 
vortex can decay by emitting both collective and single- 
particle excitations. The motion of the vortex towards 
the wall requires the presence of dissipation, as in the 
frictionless approach the velocity of the vortex core coin- 
cides with the velocity (Q) which does not contain a radial 
component. Thus, just the presence of dissipative pro- 
cesses provides a decay of the vortex state (see [^ and 
related discussion [Q of the stability of a kink state). 

The dissipation originates from the scattering of ele- 
mentary excitations by the vortex and is related to the 
friction force F in Eq.(^), which is nothing else than the 
momentum transferred from the excitations to the vor- 
tex per unit time. This force can be decomposed into 
longitudinal and transverse components: 



-Dvi — D'u X k/k, 



(5) 



where u = v^ — v„, v„ is the velocity of the normal 
component, and D, D' are longitudinal and transverse 
friction coefficients, respectively. In a static trap u = v^, 
as the normal component is at rest (v„ = 0). The friction 
force has been investigated in relation to the attenuation 
of the second sound in superfluid "^He, where the most 
important is the transverse component p3[-p5[ (see also 
|19[ for review). For a straight infinite vortex line (par- 
allel to the z-axis), a general expression for the friction 
force in a homogeneous superfluid is obtained in terms of 
the scattering amplitude /(k, k') |2^]: 
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h{ku) / (k - k' 
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6{k, - ki) 



l/(k,k')|^ 



d^k d^k' 



(27r)3 (27r)3 



- [u X k]p„ 



(6) 



Here p„ is the local mass density of the normal com- 
ponent, k, k' are the wave vectors of the incident and 
scattered excitations, n{Ek) = {exp{Ek/T) — l)~^are the 
Bose occupation numbers for the excitations, Ek is the 
excitation energy, and T the gas teniperature. Compar- 
ing the second terms of Eqs. (^ and (|^), one immediately 
arrives at the universal expression for the transverse fric- 
tion coefficient: D' = npn, assuming that the first term 
of Eq.(|) does not contribute to D' ||,|5|. 

We now turn to our analysis of the dissipative dynam- 
ics of the vortex state in a non-rotating trap, related to 
the motion of the vortex core (line) to the border of the 
condensate. This motion occurs on top of small oscilla- 
tions of the vortex filament and a slow drift (Q) of the 
vortex core. The radial component of the velocity of 
the vortex core is determined by the longitudinal friction 
coefficient D. For finding these quantities in dilute Bose- 
condensed gases, the analysis of [p3^pq,19| can only be 
used at very low temperatures (T <^ iJl)[ where the num- 
ber of thermal excitations is very small and, hence, the 
longitudinal friction force is extremely weak. 



The situation is drastically different in the temperature 
range T ^ fj., which is the most interesting for trapped 
Bose-condensed gases. We will consider the limit T ^ fi 
and first analyze how the vortex scatters excitations with 
energies Ek ^ fi. These excitations are single particles, 
and their De Broglie wave length is much smaller than 
the spatial size R of the condensate. The most important 
is the interaction of the excitations with the vortex at 
distances from the vortex line r ^ a <^ R. Therefore, the 
corresponding friction force in a trapped condensate can 
be found in the local density approximation: We may use 
Eq.(g), derived for a homogeneous superfluid, and then 
replace the condensate density no by the Thomas-Fermi 
density profile of the trapped condensate. 

The Hamiltonian of the single-particle excitations is 
h ii^ /2m, + 2nQ{r)g — ^, where the second term originates 
from the mean- filed interparticle interaction, with no(r) 
is the density of the vortex state, g = Anfi Qsc/m, Qsc is 
the scattering length, and n = no{oo)g (no(oo) = uq). 
For r ^ oo we have i? (k, r) — h k^/2m -I- /i. Hence, the 
interaction Hamiltonian responsible for the scattering of 
excitations from the vortex can be written as 

Hint = 2 [no{r)g - n] . 

For the vortex charge Z = 1, at distances r 't^ a the 
interaction Hamiltonian Hint ~ —2fj,. For r ^ a we 
have no(r) w (^ — Ti^ /2mr'^)/g, and Hint ~ —h^/mr^. 
The scattering amplitude in Eq.(|6|) can be written as 
/(k,k') = 2TTS{k^ ~ fc;)/(k,k'), where the 2D scatter- 
ing amplitude in the Born approximation is given by 



/(k,k')= d\H,nt{r)e 



?qr 



Here q = k— k' is the momentum transferred from the ex- 
citation to the vortex. As the amplitude / only depends 
on |q|, the first term in Eq.ffl) is purely longitudinal. 

For qa ^ 1, which corresponds to small angle scatter- 
ing, from Eq.(R) we obtain / '^ (?i^/m) log(l/ga). For 
qa > 1 we find \f{q)\^ r^( h^ /m)"^ sin^ {qa - 7r/4)/(ag)3. 
Using these results in Eq. (p) , we see that the main contri- 
bution to the integral over momenta comes from energies 
Ek satisfying the inequality fi ^ Ek <^ T. A direct cal- 
culation of the longitudinal friction coefficient gives 



D « Kpn{T){nog/T)'/^ 
where the density of the normal component 



(8) 



Pn 



dn n 
-P 



(Pp 



dEp"^ {2TThf 



0.1- 



^5/2y3/2 



A collective character of excitations with energies Ek ^ /i 
can influence the numerical coefficient in Eq.(g), and for 
this reason we did not present the exact value of this 
coefficient in the single-particle approximation. 

The coefficient D (x T, and Eq. (p[) can be rewritten 
as -D oc huQ^, where the quantity 






(9) 



is a small parameter of the finite-temperature perturba- 
tion theory at T ^ /i. The inequality ^ ^ 1 remains 
valid even near the BEC transition temperature, except 
the region of critical fluctuations [p6| . 

Relying on Eq. (^) for the longitudinal friction force, we 
consider the motion of the vortex line to the border of the 
condensate in a static trap, where the normal component 
is at rest. Assuming a small friction in Eqs.(P and (||), 
for finding a friction-induced small quantity v^ — V5 we 
only retain the terms linear in the dissipation coefficients 
D and D' . Then we obtain the equation 

Ps[{-VL - vs) X /t] = -Dvs - D'[vs X k]/k 



Vl ^ ^ ^L 



Vr \K, X 



which has a solution of the form vl 
r]/ nr. For the radial [v'^^^) and tangential (w^"^-*) compo- 
nents of the velocity of the vortex line we find 



,{r) 



„(</-) 



Dvs/psK, 

vs{1-D'/p,k). 



(10) 



From Eqs.([lO|) it is clear that the radial motion of the 
vortex is governed by the value of the longitudinal fric- 
tion coefficient, whereas the transverse friction (lordan- 
skii force) simply slows down the drift velocity M) of the 

(r) 

vortex. The radial velocity v)^ «C ws, which is guaran- 
teed by the inequality (||). 

The time dependence of the distance xq of the vortex 
line from the axis of a cylindrical trap follows from the 



(7) equation of radial motion for the vortex, dxo/dt 



,('■) 



With Eq.(p^ for v[ and Eq.(^ for vs, for the charac- 
teristic time of motion of the vortex from the center of 
the trap to the border we obtain 



dxQm{R^ - xDps /noffy/2 

flXQPn \ T ) 



(11) 



where Xmin is the initial displacement of the vortex line 
from the axis of the trap. The vortex velocity is the 
smallest near the trap center, and the main contribution 
to the integral in Eq.(pT]) comes from distances xq ^ R. 
Therefore, with logarithmic accuracy, we can neglect a;o 
in the numerator of the integrand and put ps = /Os(0), 
no = no(0). Then, Eq.(|l|) yields 



-h^ y—) i^^(^/^™-). 



(12) 



and 



This result is insensitive to the exact value of a;„ 
we can put Xmin ~ a. 

Once the vortex reaches the border of the condensate, 
it immediately decays. Hence, the time r can be regarded 
as a characteristic life-time of the vortex state in a static 
trap. Interestingly, the decay rate can be written as 



^0 , 3 \l/2 /^^ 



(13) 



where riom is the maximum condensate density, and Eq ~ 
h /mB? is the energy of excitation corresponding to the 
motion of the vortex core with respect to the rest of the 
condensate (excitation with negative energy, found in the 
recent calculations |2^j2|,|l^,|l|,0 ) . Eq.(^ is similar to 
the damping rate of low-energy excitations of a trapped 
condensate, found beyond the mean-field approach 
Both rates are proportional to the small parameter ^ 

For Rb and Na condensates at densities tiq ^ l{y^^cm 
and temperatures 100 ^ T < 500 nK, in the static traps 
with frequencies 10 ;S w < 100 Hz the life-time r of the 
vortex state ranges from 0.1 to 10 s. This range of times is 
relevant for experimental studies of the dissipative vortex 
dynamics. 

A proposed way of identifying the presence of a vortex 
state in a trapped Bose-condensed gas assumes switching 
off the trap and observing a ballistically expanding gas 
sample [O]. As follows from the numerical simulations 
ijlll , at zero temperature the expansion of a condensate 
with a vortex occurs along the lines of the scaling theory 
p9| , p0[ . The shape of the Bose-condensed state is nearly 
preserved and its spatial size is increasing. Due to ex- 
pansion the density of the condensate decreases, and the 
size of the vortex core increases to match the instanta- 
neous value of the healing length. This should allow one 
to detect the vortex through the observation of a hole in 
the density profile of the condensate. 

It is important to emphasize that at temperatures 
T ^ II the thermal cloud will expand with the ther- 
mal velocity vt ^ \/T jra which is much larger than 
the expansion velocity of the condensate (the latter is 
of order the sound velocity Cg). Therefore, after a short 
time RjvT the thermal component files away, and the 
dissipation-induced motion of the vortex core ceases. Ac- 
cordingly, the expansion of the Bose-condensed state will 
be essentially the same as that at zero temperature. This 
means that the relative displacement of the vortex core 
from the trap center practically remains the same as be- 
fore switching off the trap. Therefore, the dissipative 
dynamics of the vortex state in the initial static trap, i.e. 
the motion of the vortex core towards the border of the 
condensate, can be studied by switching off the trap at 
different times and visualizing the position of the vortex 
core in a ballistically expanding condensate. 

In conclusion, we have developed a theory of dissipative 
dynamics of a vortex state in a trapped Bose-condensed 
gas at finite temperatures and calculated the decay time 
of the vortex with charge equal to 1 in a static trap. Our 
theory can be further developed to analyze the motion of 
vortices in rotating traps and, in particular, to calculate 
a characteristic time of the formation of the vortex state 
in a trap rotating with supercritical frequency. 
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